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摘要 ^ AURIC 是 由 美国 计算 物理 公司 CPI 与 空军 Phillips 实验 室 联 合 开 发 的 中 高 层 大 气 紫 “| 关键 词 
外 -可 见 光 - 近 红外 气 辉 辐射 传输 模型 , 是 目前 唯一 用 来 进行 中 高 层 大 气 气 辉 辐 笑 模 拟 研究 的 | 人 全 ep 
通用 模型 . 基于 MODTRAN 模型 的 理论 , AURIC 可 以 进行 80 km 高 度 以 上 的 辐射 传输 模拟 并 Non 
将 辐射 波段 扩展 到 远 紫 外 波段 (80 nm). E 81, CPI 公司 只 面向 全 球 发 布 了 AURIC v2 软件 封 | 体 辐 射 率 


装 包 ,其 只 能 进行 2000 年 以 前 的 单 点 计算 , 不 适用 于 批量 的 中 高 层 大 气 辐射 模拟 ,更 不 能 用 
于 星 上 大 气 成 分 的 批量 反 演 . 本 文 利用 Matlab 对 AURIC v1.2 地 磁 参 数 模块 、 大 气 电子 密度 
模块 、 大 气温 度 及 各 种 中 性 气体 成 分 密度 模块 进行 了 替代 ,将 更 新 模块 与 原 有 的 辐射 计算 模 
块 相 结 合 , 将 AURIC v1.2 更 新 为 AURIC-2012 模型 ， 其 可 以 批量 地 进行 全 球 的 大 气 辐 射 传输 
模拟 ,可 以 与 星 载 测量 数据 相 结 合 进 行 中 高 层 大 气 成 分 的 批量 反 演 ， 如 O/N> .电子 密度 等 ， 同 
时 也 为 模型 辐射 计算 模块 的 进一步 改进 和 辐射 机 制 的 参数 更 新 奠定 基础 . 基于 AURIC-2012 
模型 ,进行 了 气 辉 临 边 柱 辐 射 强 度 模 拟 和 体 辐射 率 计算 并 将 结果 分 别 与 GUVI 柱 辐 射 强度 
和 TIDI 体 辐 射 率 实测 值 进 行 了 比较 , 得 到 两 者 峰值 的 模拟 平均 相对 误差 都 小 于 20%. 最 后 ， 
利用 AURIC-2012 对 气 辉 临 这 柱 辐 射 强度 随 纬 度 和 高 度 的 分 布 进 行 了 二 维 模拟 , 并 基于 模拟 
对 上 怪 、 夜 气 辉 辐射 强度 分 布 特性 及 影响 因素 进行 了 分 析 . 


气 辉 是 中 高 层 大 气 一 种 重要 的 自然 发光 现象 ， ”外 气 辉 . 在 各 种 谱 线 气 涟 中 , 不 同 的 气 辉 分 布 高 度 不 


其 主要 是 由 太阳 紫外 辐射 直接 或 间接 激发 大 气 分 子 
或 原子 而 产生 的 . 气 辉 主要 分 布 在 80 km 以 上 高 度 ， 
mS MAARE, 并 且 剧 夜 一 直 存 在. 按照 发 生 
KETE, SJEN ADNE AIE, KAEM E A OE; 
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同 ， 辐 射 强度 也 不 同 ; qf E — X E T vm SE BG T 
理 位 置 和 时 间 也 在 不 断 的 变化 ,强度 的 变化 与 大 气 
的 化 学 过 程 和 动力 学 过 程 存 在 着 密切 的 关系 ,这 表 
明 气 辉 是 反应 这 些 物 理 过 程 和 物理 参数 的 媒介 ， 因 
而 被 蕉 为 中 高 层 大 气 的 “指纹 ” 因此 , 利用 模型 对 气 
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辉 辐射 的 模拟 计算 ， 有 助 于 深入 了 解 气 辉 的 辐 冉 特 
lE, 并 对 了 解 电 离 层 的 化 学 、 动 力学 及 热力 学 等 过 程 
提供 参考 数据 . 同时 , 利用 模型 与 气 辉 的 星 载 测量 值 
相 结 合 ,还 可 以 对 电离 层 大 气 成 分 进行 及 演 (如 O/N;， 
电子 密度 等 )， 为 科学 研究 和 军事 活动 等 提供 观测 
数据 . 


HU. ES bn EIOS HP RACK] URS R A 
一 般 只 基于 一 种 或 者 两 种 气 辉 进行 它们 的 体 辐 射 率 
或 强度 的 模拟 (Solomon 和 Abreu，1989; Richards 和 
Torre, 1988; Richards 等 , 1994; Murtagh 和 Donal, 1989; 
Murtagh 等 ,1990; Witasse 等 , 1999; Estes, 2000; L 
等 , 2014)， 而 综合 起 来 进行 中 高 层 大 气 气 辉 辐 射 传输 
模拟 的 通用 模型 只 有 大 气 紫 外 辐射 模型 AURIC 
(Atmospheric Ultraviolet Radiance Integrated Code). 
该 模型 是 由 美国 计算 物理 公司 CPI (Computational 
Physics, IncJ 和 空军 Phillips 实 验 室 经 过 十 多 年 的 研究 
于 2002 年 联合 发 布 的 ， 以 用 于 中 高 层 大 气 辐射 传输 
模拟 (Link 等 ，1992;，Strickland 等 ，1999，Evans 等 ， 
2002). 1999 年 , Strickland 等 第 一 次 利用 AURIC 模 型 进 
行 气 辉 辐射 强度 、 体 辐射 等 的 计算 模拟 ， 并 将 模拟 结 
条 与 火箭 、 卫 星 等 实测 数据 进行 了 对 比分 析 . 之 后 ， 
AURIC 模 型 逐渐 的 被 应 用 于 大 气 成 分 的 反 演 等 方面 . 

Bishop 等 (2003) 对 AURIC v1.0 极 紫外 波段 范 
进行 了 进一步 扩展 ， 并 基于 该 模型 对 霍 普 金 斯 紫外 
望远镜 (HUT) 观 测 的 极 紫 外 数据 进行 了 分 析 ; Stephan 
等 (2004) 利 用 AURIC 对 赤道 地 区 热 层 人 磁 暴 时 紫外 辐 
射 情况 进行 了 研究 ; Strickland 等 (2004) 利 用 AURIC 建 
立 的 查找 表 得 到 了 O/Ns 与 135.6/LBHs 的 关系 ,并 结 
合 GUVI 数 据 进 行 全 球 热 层 O/N, 的 季节 变化 分 析 ; 
Bobik 等 (2013) 利 用 AURIC 对 300~400 nm 的 夜 气 辉 辐 
射 率 进 行 计 算 ， 并 用 于 南大 西洋 异常 区 的 紫外 辐射 
背景 研究 ; Evans 等 (2013) 结 合 火 星 大 气 的 气体 成 分 
及 数据 对 AURIC 模 型 进行 了 和 更改， 以 用 于 对 火星 大 
气 辐射 的 研究 . 在 国内 , 近 几 年 学 者 基于 AURIC 也 进 
行 了 一 些 研究 ， 绢 圣 峰 等 (2012) 利 用 AURIC 模 型 对 
135.6 nm 和 LBH 昼 气 辉 进行 了 模拟 研究 ， 并 将 其 用 于 
热 层 O/N 计 算 ; 张 永超 等 (2014) 基 于 AURIC 模 型 进行 
了 LBH 计 算 方法 的 改进 ， 引 入 了 观测 方位 角 ， 并 考虑 
了 大 视 场 下 太阳 天 顶 角 变化 对 结果 的 影响 . 

目前 ， 由 于 很 难 获 得 AURIC v1.2 的 源 代 码 ， 绝 
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大 部 分 研究 一 般 仍 基于 CPI 发 布 的 AURIC v1.2 软 件 
包 进 行 的 . 而 AURIC v1.2 只 是 一 个 软件 封装 包 ， 这 
使 得 其 只 能 进行 单 点 计算 ， 不 能 进行 批量 模拟 ， 且 不 
能 用 于 星 载 大 气 成 分 的 批量 数据 反 演 . 此 外 ，AURIC 
v1.2 电 子 密度 计算 是 基于 FAIM 模 型 ,而 中 性 气体 密 
度 和 温度 计算 基于 NRLMSISE-00 和 SHARC 模 型 都 
是 基于 2000 年 及 以 前 的 数据 库 . 因此 ， 本 文 利用 国际 
参考 电离 层 模型 IRI-2012 电 子 密度 计算 模块 和 
DTM-2013 大 气温 度 、 中 性 气体 密度 计算 模块 对 
AURIC v1.2 进 行 了 模块 更 新 ， 并 基于 Matlab 将 其 集 
成 为 AURIC-2012 程 序 包 ,可 循环 进行 全 球 任 意 时 间 
的 气 辉 辐射 模拟 , 并 可 用 于 中 高 层 大 气 成 分 的 星 载 
反 演 . 


1 原理 与 方法 


在 AURIC v1.2 中 ， 昼 气 辉 计算 模型 主要 有 14 个 
模块 ， 分 别 为 地 磁 参 数 GEOPARM、 大 气 模式 模块 
ATMOS、 电 离 层 模式 模块 IONOS、 太 阳光 谱 模 块 
SOLAR、 倾 角 柱 密度 计算 模块 COLDEN、 光 电子 源 
模块 PESOURCE、 光 电子 通 量 模 块 PEFLUX、 电 子 倍 
撞 模块 E-IMPACT、 化 学 反应 模块 DAYCHEM、 多 元 
汇合 模块 MERGEVER、 观 测 视线 柱 密度 计算 模块 
LOSDEN、 和 辐射 传输 模块 RADTRANS、 柱 辐射 率 计 
算 模块 LOSINT 和 柱 辐 射 率 结果 存储 模块 
MERGEINT. 夜 气 辉 计算 模型 则 主要 有 8 个 模块 ,分 
别 为 地 磁 参 数 GEOPARM、 大 气 模式 模块 ATMOS、 
电离 层 模式 模块 IONOS、 夜 气 辉 辐射 计算 模块 
NITEGLO、 观 测 视线 柱 密度 计算 模块 LOSDEN 、 辑 
射 传输 模块 RADTRANS、 柱 辐射 率 计 算 模块 LOSINT 
和 柱 辐 射 率 结果 存储 模块 MERGEINT. 这 些 模块 封 
装 在 同一 个 bat 执 行文 件 中 ,每 个 模块 承担 一 种 计算 
功能 ， 模块 之 间 通 过 参数 或 者 文件 相互 连接 ， 因 此 
AURIC 可 以 依次 调用 不 同 的 计算 模块 来 分 别 进行 入 
气 辉 和 夜 气 辉 的 模拟 计算 . 

由 于 AURIC v1.2 只 能 进行 单 点 的 气 辉 模拟 ,不 
适用 于 大 量 数 据 的 仿真 模拟 及 星 上 大 气 成 分 反 演 . 
因此 ， 本 文 对 地 磁 参 数 GEOPARM 模 块 、 电 离 层 
IONOS 模 块 及 大 气 ATMOS 模 块 进 行 了 更 新 奉 代 ， 改 
进 内 容 及 步骤 如 图 1 所 示 . 
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磁 经 纬度 \ 磁 倾角 


太阳 天 顶 角 
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改进 2 [温度 及 中 性 密度 — Conn ) 


改进 1 


as [a CR 


图 1 AURIC-2012 改进 内 容 及 步骤 


图 1 中 ,改进 1 为 AURIC v1.2 中 GEOPARM 模 块 
参数 的 更 新 替代 ,主要 包括 地 磁 经 纬度 、 磁 倾角 、 太 
HAMA, 当地 时 间 、 太 阳 指 数 Fio7 及 地 磁 指 数 Ap. 

这 些 参数 只 更 新 到 1999 年 12 月 31 日 , 这 使 得 AURIC 
不 适用 于 目前 观测 数据 的 应 用 . 因此 本 次 更 新 对 这 
些 参数 进行 了 重新 计算 和 模块 代替 . 其 中 ， 磁 倾角 是 
由 国际 地 磁 参 考 模型 (Finlay 等 , 2010) 计 算得 到 的 ; 大 
虽 天 顶 角 的 计算 是 Msn a n 
置 算法 ; 太阳 Fio7 指 数 和 地 磁 Ap 指 数 的 计算 则 是 
于 NOAA 提 供 的 实测 数据 库 EE 
gov/STP/GEOMAGNETIC_DATA/INDICES/KP_AP). 

改进 2 为 大 气温 度 和 中 性 气体 密度 廓 线 的 计算 . 
AURIC v1.2 中 ATMOS 模 块 的 计算 是 基于 2000 年 以 前 
的 数据 库 进 行 的 ， 这 不 适用 于 目前 的 气 辉 模 拟 和 大 
气 成 分 星 上 反 演 . 因此 本 次 更 新 基于 DTM-2013 模 型 
进行 大 气温 度 和 中 性 气体 密度 的 计算 ， 并 将 其 作为 
AURIC-2012 中 ATMOS 的 更 新 模块 . DTM 模 型 是 基于 
卫星 实测 数据 的 半 经 验 模 型 ,主要 用 于 大 气温 度 和 
中 性 气体 密度 的 模拟 计算 ， 其 中 中 性 气体 主要 包括 
H, He, O, 0; 和 Ns. 目前, DTM 模 型 先后 经 历 了 6 个 升 
级 版 本 ,其 是 由 Barlier 等 (1978) 利 用 从 卫星 上 忠 力 数据 
反 演 的 热 层 大 气温 度 和 密度 数据 首次 开发 的 ， 之 后 
经 过 DTM-1994, DTM-2000, DTM-2009, DTM-2012 
及 DTM-2013 版 本 (Berger 等 , 1998; Bruinsma 等 , 2003, 
2004, 2012). 因此 ， 本 文 利 用 其 最 新 的 版 本 DTM- 
20134& Ë AURIC v1.2 的 ATMOS 模 块 进行 温度 和 中 性 
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气体 密度 的 计算 ， 从 而 实现 ATMOS 模 块 的 更 新 . 

改进 3 为 电子 密度 廓 线 的 计算 . AURIC v1.2 中 
IONOS 模 块 的 计算 模拟 是 基于 1989 年 以 前 的 实测 电 
子 密度 数据 库 进 行 的 ， 这 不 适用 于 现在 的 气 辉 模拟 
和 大 气 成 分 星 上 反 演 . 因此 本 次 更 新 是 利用 国际 参 
考 电离 层 IRI2012 代 蔡 原 IJONOS 模 块 进行 电子 密度 
的 计算 ， 并 将 其 作为 AURIC-2012 的 IONOS 模 块 . IRI 
模型 是 在 20 世 纪 60 年 代 由 空间 研究 委员 会 (COSPAR) 
和 国际 无 线 电 科 学 联盟 (URSD 共 同 发 起 建立 的 电离 
层 模型 ， 其 是 基于 世界 范围 内 实测 数据 的 经 验 预报 
模型 . 它 综合 了 全 球 180 多 个 电离 层 探测 站 的 资料 ， 
主要 描述 了 海拔 高 度 50~1500 km 范围 内 电离 层 的 电 
子 密度 、 电 子 温度 、 离 子 成 分 及 离子 温度 等 重要 参数 
(Bilitza 和 Reinisch, 2008; Bilitza 等 , 2014). 因此 ， 本 文 
利用 其 最 新 的 版 本 了 IRL2012 代 蔡 AURIC v1.2 中 
IONOS 模块 进行 电子 密度 廓 线 计 算 ， 从 而 实现 
IONOS 模 块 的 更 新 . 

改进 后 的 AURIC-2012 模 型 主要 输入 参数 有 三 
个 ,分 别 是 地 理 纬 度 、 地 理 经 度 和 世界 时 (UTC)， 而 
AURIC-2012 模 型 的 输出 参数 主要 为 各 种 气 辉 的 柱 辐 


射 强度 、 体 辐射 率 等 参数 . AURIC-2012 模 型 计算 流程 
如 图 2 所 示 . 
输入 参数 


图 2 AURIC-2012 计算 模拟 流程 图 


2 模型 验证 与 分 析 


基于 AURIC-2012， 本 文 对 GUVI 气 辉 观测 进行 
了 仿真 模拟 , 包括 130.4，135.6 nm 和 LBH 波 段 气 辉 ， 
并 将 模拟 结果 与 实测 数据 进行 了 比较 . 以 2002 年 3 月 
20 日 地 理 坐 标 为 (31.78°"N，138.4*W) 的 GUVI 数 据 为 
例 ， 对 比 结果 如 图 3 所 示 . 图 中 ，130.4，135.6 nm 和 
LBHJE Ez PE 48 4] vti BE B UL ELE; GU V ELSE WE BE ne BE 
的 变化 规律 ( 纵 坐 标 ) 较 为 一 致 ， 同 一 位 置 同 一 高 度 处 


AURIC-2012 


(a) 


度 (km) 


5 
[e 


一 
n] 


高 度 (km) 


高 度 (km) 


度 (km) 
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两 者 的 辐射 强度 值 相近 ， 其 中 3 个 波段 气 辉 在 峰值 区 
域 的 模拟 相对 误差 分 别 为 21.93%(130.4 nm), 
22.57%(135.6 nm) 和 11.32%(LBH), 平均 相对 误差 为 
18.6%， 且 峰值 模拟 值 与 GUVI 峰 值 实测 值 处 于 同 


将 AURIC-2012 计 算得 到 的 0,(0-0) 波 段 气 辉 体 辐 
fj 3 5 TIDIDUL JU Zi dg xtír fo, un EAT. 
AURIC-2012 体 辐射 率 计 算 模拟 值 随 高 度 分 布 规律 与 
TIDI 观 测 数 据 相 吻 合 ，4 次 模拟 峰值 区 域 (3 个 采样 


GUVI (x104) 
4 


图 3 130.4 nm(a), 135.6 nm(b) 和 LBH 波段 (@) 气 辉 辐 射 强度 对 比 (单位 : 瑞 利 ) 
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图 4 体 辐射 率 模拟 值 与 TIDI 实测 值 比较 结果 
(a) 2005 年 9 月 2 日 ;(b)2008 年 5 月 17 日 ;(c)2009 年 5 月 17 日 ;(d)2010 年 1 月 5 日 


点 ) 相 对 误差 分 别 为 8.56%( 图 4(a))，12.49%( 图 4(b))， 
24.67%( 图 4(c)) 和 16.06%( 图 4(d)), 平均 相对 误差 为 
15.45%， 旦 峰值 高 度 模 拟 值 与 TIDI 实 测 值 非常 一 致 . 


3 气 辉 辐射 模拟 分 析 
3.1 气 辉 随 纬 度 和 高 度 变化 模拟 

基于 AURIC 的 更 新 改进 ,本 文 利用 AURIC-2012 
进行 了 随 纬度 和 高 度 变化 的 临 边 柱 辐 射 强度 模拟 ， 
模拟 日 期 选择 多 个 年 份 的 春分 点 附近 ， 以 116?E( 跨 北 
AO by ti, 纬度 窗 盖 范围 为 90*S~90°N, 模拟 
高 度 为 80~600 km. 由 于 AURIC-2012 只 进行 了 5 种 主 
要 中 性 气体 密度 (HE，He，O，O 和 N?) 的 计算 更 新 ， 
此 本 文 的 气 辉 模拟 分 析 主 要 包括 氧 原子 气 辉 谱 线 
130.4, 135.6, 557.7, 630.0, 636.4, 777.4 和 844.6 nm. F 
A 4 Y ^U TÉ ZLBHCE ^U) UL SUL TA 
谱 线 O,(0-0)( 仅 夜 气 辉 )， 这些 气 辉 也 是 用 于 电离 层 探 
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测 和 研究 最 主要 的 谱 线 气 辉 . 模拟 结果 以 2013 年 3 月 
20 日 为 例 (Fi07=106.7, Aps=12, Apx=9)， 如 图 5~6 所 示 
(图 5 为 生气 辉 ， 图 6 为 夜 气 辉 ). 

图 5 为 生气 辉 柱 辐射 强度 的 模拟 结果 ,分别 为 
130.4, 135.6, 630.0, 636.4, 777.4, 844.6, 557.7 nm 以 及 
LBH 波 段 气 辉 . 由 图 5 可 知 , EU dS EI JLSE TS os PT 
有 纬度 ， 其 中 在 赤道 附近 最 强 ， 随 着 纬度 的 增加 气 瘤 
辐射 强度 逐渐 减弱 . 而 昼 气 辉 辐 射 随 高 度 一 般 先 增 
大 后 减 小 ,但 不 同 谱 线 气 辉 的 主要 分 布 高 度 和 变化 
快慢 规律 不 同 . 130.4，135.6，630.0，636.4，777.4 和 
844.6 nm 气 辉 辐 射 分 布 高 上 度 相 对 较 高 ， 随 高 上 度 变 化 相 
对 缓慢 ; 而 557.7 nm 和 LBH 气 辉 的 分 布 高 度 相 对 较 
低 ， 随 高 度 的 变化 较 快 . 每 个 谱 线 生气 辉 具 体 的 分 布 
和 变化 规律 如 下 : 

(1) 130.4 nm 昼 气 瘤 较 强 辐射 主要 分 布 在 300 km 
高 度 以 下 ， 辐 射 强度 最 大 可 达到 几 万 瑞 利 ; 当 >300 
km 时 ， 随 着 高 度 的 增加 ,辐射 强度 逐渐 减 小 到 几 干 
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图 5 昼 气 辉 临 边 柱 辐射 强度 模拟 结果 (单位 : 瑞 利 ) 


瑞 利 (<600 km). 

(2) 135.6 nm 有 昼 气 辉 辐射 也 主要 分 布 在 300 km 以 
下 ,辐射 强 度 最 大 为 儿 王 瑞 利 ; 在 >300 km 时 , 气 辉 
辐射 迅速 减 小 到 几 百 瑞 利 ; 而 当 >500 kmit, 气 辉 则 
变 得 非常 微弱 ， 仅 为 几 个 瑞 利 . 

(3) 630.0 和 636.4 nm 气 辉 的 辐射 机 理 一 致 ， 两 者 
的 强度 变化 规律 完全 相同 , 但 由 于 跃迁 能 级 不 同 ， 
此 两 种 气 辉 的 辐射 强度 不 同 (1uaoo16364=3.09).， 其 中 ， 
630.0 nm& ^UE EI 4E 1E300 km 以 下 ， 辐 射 强 
度 达 到 几 万 瑞 利 ; 342300 km 时 ， 随 着 高 度 的 增加 ， 
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辐射 强度 逐渐 减 小 ; 15252500 kmit, 气 辉 则 变 得 非 
Tí dus, DUNSJVT SI. 

(4) 777.4 nm E ^UE AH] 35 927 4h 16250 km 以 下 ， 
辐射 强度 为 几 千 瑞 利 ; 当 >250 km 时 ， 随 着 高 度 的 增 
加 气 辉 辐射 强度 逐渐 减 小 ; 而 当 >400 km 时 , SIE 
射 强 度 只 为 几 个 瑞 利 . 

(5) 844.6 nm 昼 气 辉 辐 射 变化 规律 与 777.4 nm 一 
致 ， 但 是 辐射 强度 相对 较 强 ,这 是 由 于 844.6 nna A 
辉 除了 由 氧 原子 与 光电 子 复合 产生 之 外 ,还 来 自 于 
氧气 分 子 及 氧 原 子 与 电子 的 反应 贡献 ( 表 1). 
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(6) 557.7 nm 气 辉 辐射 主要 分 布 在 200 km 以 下 ， 
其 辐射 强度 较 大 ,达到 几 十 万 瑞 利 ; 当 >200 km 时 ， 
气 溜 辐射 强度 迅速 减 小 到 几 和 王 瑞 利 ; 而 当 >400 km 
时 ， 气 瘤 辐射 只 为 几 个 瑞 利 . 

(7) LBH 气 辉 与 557.7 nm 气 辉 辐射 分 布 在 同一 
RE, 8T EJL H Fi A. 不 同 于 其 他 谱 线 气 辉 ， 
LBH 和 气 瘤 为 远 紫 外 宽 谱 段 气 辉 ,波段 覆盖 范围 为 
140-190 nm. 

图 6 为 夜 气 辉 临 边 柱 辐射 强度 的 模拟 结果 ,分 别 
为 130.4, 135.6, 630.0, 636.4, 777.4, 844.6, 557.7 nm 以 


基于 紫外 辐射 传输 模型 AURIC-2012 的 气 辉 辐射 模拟 
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及 0;(0-0) 波 段 气 辉 . FH AURICAP BE 3E £1 88 3-008 
(90°< 太 阳 天 顶 角 <110°) 计 算 ， 因此 3 月 20 日 夜 气 辉 辐 
射 强 度 计 算 主要 位 于 纬度 (70"S~70*N) 范 围 内 . 由 图 
5~6 可 知 ， 夜 气 辉 的 水 平 变 化 规律 与 骨气 辉 一 致 ， 在 
赤道 附近 较 强 ， 且 随 着 纬度 的 增加 气 辉 辐射 强度 志 
速 减弱 . BARTET, ROC OPERI RER, Zh 
度 分 布 范围 也 相对 较 小 ， 除 557.7 nm 和 OO;(0-0) 波 段 气 
辉 外 ， 其 他 夜 气 辉 辐射 主要 分 布 在 中 低 纬 上 度 地 区 . 此 
外 , 夜 气 辉 随 高 度 的 分 布 和 变化 规律 与 生气 辉 也 不 
一 样 , 130.4, 135.6, 630.0, 636.4, 777.4 和 844.6 nm 夜 气 
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图 6 夜 气 辉 临 边 柱 辐射 强度 模拟 结果 (单位 : 瑞 利 ) 
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辉 的 辐射 覆盖 范围 较 广 ,变化 相对 缓慢 ， 且 峰值 高 度 
明显 较 高 ; 而 557.7 nm 和 O5(0-0) 夜 气 辉 分 布 的 峰值 高 
度 相 对 于 昼 气 辉 则 变 得 更 低 (峰值 为 ~95 km), 辐射 高 
度 范 围 也 更 窗 ,， 仅 为 ~40 km. 每 个 谱 线 夜 气 辉 具体 的 
分 布 和 变化 规律 如 下 : 

(1) 130.4 nm 夜 气 辉 辐射 相对 于 昼 气 辉 要 弱 很 多 ， 
最 大 只 有 几 百 瑞 利 ， 主 要 分 布 在 低 纬 度 地 区 500 km 
高 度 以 下 ， 中 高 纬度 气 辉 辐射 强度 只 有 几 个 瑞 利 . 

(2) 135.6 nm 夜 气 辉 辐射 主要 分 布 在 低 纬度 
100-500 km 高 度 范围 内 ， 最 大 峰值 出 现在 ~300 km 高 
度 ， 辐 射 强度 为 几 百 瑞 利 ， 而 当 >500 kmi, AXE M 
变 得 非常 微弱 . 

(3) 630.0 和 636.4 nm ÀX ^ XE 4m S] (Iaso.o/ Las. = 
3.09) 在 中 低 纬度 ~100~400 km 高 度 范围 内 辐射 较 强 ， 
其 中 630.0 nm 气 辉 辐射 达 几 千 瑞 利 ; 在 高 纬 地 区 依然 
4 ^U. ZJJLE SRI; 无 论 中 低 纬 还 是 高 纬 地 
区 ， 最 大 峰值 一 般 出 现在 ~250 km 高 度 ， 而 当 >400 
km 时 ， 气 辉 则 变 得 非常 微弱 . 

(4) 777.4 和 844.6 nm 夜 气 辉 辐射 主要 分 布 在 中 低 
纬 ~150~450 km 高 度 范 围 内 ， 最 大 峰值 出 现在 
300-350 km 高 度 ， 辐 射 强 度 为 几 百 瑞 利 ， 而 在 高 纬 地 
区 气 辉 则 变 得 非常 微弱 图 中 , 135.6 nm 气 辉 与 777.4 
和 844.6 nm 气 辉 辐射 变化 规律 一 怪 , 这 是 因为 3 种 夜 
气 辉 均 主要 由 光电 子 与 氧 原子 碰撞 及 离子 中 和 反应 
产生 , 但 是 跃迁 能 级 、 跃 迁 概 率 和 复合 系数 不 同 ， 所 
以 3 种 气 辉 的 辐射 强度 会 有 一 定 差异 . 

(5) 557.7 nm 夜 气 辉 辐射 覆盖 纬度 范围 较 广 ， 其 
中 在 中 低 纬 较 大 , 在 高 纬 地 区 有 所 减 小 ,但 辐射 仍然 
较 强 ， 为 几 王 瑞 利 ; [HR EI om BE en REY HUE, 位 
于 80~120 km， 强 度 最 大 可 达 几 万 瑞 利 ; 但 当 >120 
km 时 ， 气 辉 辐射 强度 迅速 减 小 并 变 得 非常 微弱 . 

(6) O02(0-0)A 波 段 夜 气 辉 辐射 分 布 规律 与 557.7 
nm 气 辉 类 似 , 但 其 辐射 强度 更 大 ， 最 大 达到 几 十 万 
瑞 利 . O2 (0-0) A 波段 气 辉 受 氧气 密度 的 影响 较 大 ， 而 
557.7 nm 气 辉 则 主要 受 氧 原子 密度 的 影响 (Strickland， 
1999), 在 ~100 km 高 度 处 氧气 密度 远 远 大 于 氧 原 子 
的 密度 ， 同 时 O;(0-0)A 波 段 气 辉 的 湾 灭 系数 远 小 于 
557.7] nm ^ X If] TE K 28 Zt (Strickland, 1999). 
JE, Ox(0-0) AJ& Ex OPE AR] Eo 1557.7 nm 气 辉 的 
辐射 . 

由 图 5 和 6 可 得 ,， 气 辉 辐射 强度 随 纬度 和 高 度 分 
布 具有 一 定 的 规律 : 
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(1) 氧 原子 、 氧 气 分子 及 毛 气 分 子 气 辉 辐射 在 赤 
道 低 纬 地 区 明显 强 于 其 他 地 区 ， 随 着 纬度 的 增加 逐 
渐 减 弱 ， 这 主要 受 电子 密度 纬度 分 布 的 影响 . 大 部 分 
气 辉 辐射 来 自 于 自由 光电 子 与 其 他 粒子 的 碰撞 或 有 反 
应 ( 表 1), 而 赤道 低 纬 地 区 电离 层 存 在 赤道 异常 ， 导 
致 该 区 域 电 子 密度 明显 高 于 其 他 高 纬 地 区 ， 因 此 气 
辉 辐射 在 该 区 域 也 明显 增强 ， 随 着 纬度 的 增加 ， 电 子 
密度 的 减 小 , 气 辉 辐 射 强度 也 随 之 减弱 . 

(2) 气 辉 辐 射 强度 随 着 高 上 度 的 增加 先 增 大 后 减 
小 . 气 辉 一 般 是 由 电子 和 其 他 至 少 一 种 粒子 进行 相 
互 作用 产生 的 ， 由 于 电子 密度 随 高 度 不 断 增 大 (<500 
km)， 而 其 他 大 部 分 粒子 随 高 度 一 般 呈 现 相 反 的 变化 
规律 (图 7)， 因此 气 辉 的 辐射 强度 随 高 度 的 增加 一 般 
先 增 大 后 减 小 ， 即 在 某 一 高 度 将 存在 一 个 辐射 最 大 值 . 

同时 由 图 5~6 可 知 ， 入 气 辉 与 夜 气 辉 的 辐射 强度 
分 布 也 存在 一 些 不 同 的 规律 : 

(1) BAHAI ERAK TAA, HARE 
盖 范 围 也 相对 较 大 . 由 于 白天 大 气 受 太阳 较 强 的 直 
接 照射 ， 电 子 的 密度 会 大 大 增加 (图 7(a))， 其 他 粒子 
也 有 不 同 程度 的 增加 (图 7(b)), 再 加 上 气 辉 的 激发 可 
以 通过 太阳 的 光 致 电离 而 直接 产生 ,， 气 辉 的 辐射 机 
制 也 变 得 相对 复杂 多 样 ， 因 此 和 昼 气 辉 的 辐射 强度 及 
分 布 范围 也 会 大 大 增加 . 

(2) 557.7 nm 夜 气 溜 辐射 高 度 范围 变 军 ， 峰 值 高 
度 相 对 较 低 . 557.7 nm 夜 气 辉 的 辐射 主要 来 自 于 氧 原 
子 的 三 体 复合 (需要 第 三 种 中 性 气体 粒子 的 参与 ， 如 
O; 和 Ns)， 因 此 其 主要 发 生 在 ~100 km 高 度 及 以 下 区 
域 . 而 557.7 nB ALAXX TARTS E T A iÉ 
撞 、 太 阳 直 接 辐 射 导 致 的 氧气 分 子 电 离 以 及 氮气 的 惰 
性 化 过 程 ， 这 些 机 制 在 100~200 km 高 度 范 围 内 均 有 
较 大 的 辐射 贡献 ， 使 100~200 km 高 度 成 为 白天 最 主 
要 的 辐射 区 域 (Upadhayaya 等 ，2006; Culot 等 ，2004; 
Thirupathaiah 等 ，2012). KE, 4X TAA, 557.7 
nm 夜 气 辉 辐射 峰值 高 度 会 变 得 较 低 ,辐射 厚 度 范 围 
变 得 较 罕 . 

(3) 130.4 nm 夜 气 辉 辐射 在 中 层 -低热 层 区 域 辐射 
强度 随 着 高 度 的 增加 变化 非常 小 , 这 是 由 氧 原 子 与 
光电 子 们 撞 及 辐射 传输 综合 作用 的 结果 ， 这 个 高 度 
范围 主要 由 电子 密度 (或 峰值 ) 的 分 布 高 度 决 定 的 
(Dymond, 2009). 

由 气 辉 的 辐射 机 制 可 知 ， 中 高 层 大 气 电 子 密度 
对 气 辉 的 辐射 强度 具有 和 较 大 的 影响 ， 是 气 辉 辐 射 分 
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表 1 
气 辉 辐射 线 


各 种 气 辉 的 辐射 机 制 ” 


辐射 机 制 
氧 原子 与 光电 子 碰撞 ( 朋 气 辉 ): O +e 2 OGS) «e 


太阳 光子 共振 散射 ( 恒 气 辉 ): O+Hv — O(3S)+hv(130.4 nm) 


OI 130.4 nm 


OI 135.6 nm 


OI 630.0 nm/636.4 nm 


OI 777.4 nm 


OI 844.6 nm 


OI 557.7 nm 


氧气 分 子 离 解 : Oe 一 O+O(3S)+e 
光 致 电离 (生气 辉 ): O?+hv 一 O+0(35) 
氧 离 子 与 电子 复合 : Of++e -> O-hv(1304 nm) 
离子 中 和 反应 : O++O- > 20+hv(130.4 nm) 
跃迁 : O(3S) 2 O(3P)+hv(130.4 nm) 


氧 离子 与 电子 复合 : Or*+e 0-hv(135.6 nm) 


氧 离子 中 和 反应 : O*- 07 — 20+hv(135.6 nm) 
氧气 分 子 离 解 : 
离 解 复合 : 0; +fe 一 O(D')+O 
JETA: e «02e «O(D') 

光 致 电离 ( 昼 气 辉 ): O, «hv 2 O(P?) - O(D') 


p. 


O,+e 一 O+O+e *hv(135.6 nm) 


离 解 复 合 : NO* +e” >0(D')+N 

禁 戒 跃迁 : 

氧 离子 与 电子 复合 : Ot -e—0-hv(7774 nm) 
离子 中 和 反应 : O01 +O- — 20 +hv(777.4 nm) 
ETETE): Ore 一 OGP)+e 
离 解 复合 ( 昼 气 辉 ): O, +e 一 OGP)+O+e- 


O (D') — O + hv(630.0 nm) 


氧 离子 与 电子 复合 : 0 +e — O +hv(844.6 nm) 


离子 中 和 反应 : 0 +O- — 20 +hv(844.6 nm) 
电子 跃迁 : O(3P) O(3S)- hv(844.6 nm) 


氧 原 子 与 电子 碰撞 : O+e 一 O(S)J+e- 
光 致 电离 (生气 辉 ): O, hv  O« O(S) 
离 解 复合 : 0,* +e 0+0) 


氮气 惰性 化 : N (AS 0*0 N, &O(IS) 


三 体 复合 :O+O+M 一 0;+M M=N,,0，0;+0 一 O(S')+0， 


O;-A (0-0) 762.0 nm 


N2 LBH 


a) 5| É Strickland 等 (1999), Ghosh(2002) 8I Witasse 等 (1999) 


布 的 最 主要 影响 因素 ,尤其 对 于 辐射 机 制 相对 简单 
的 夜 气 辉 ， 电 子 密 度 的 变化 将 直接 影响 它们 的 辐射 
变化 . 图 8 为 不 同 电 子 密度 时 135.6 和 630.0 nm 气 辉 体 
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禁 戒 跃迁 
三 体 复 合 : 0+0+M>0,(b')+M M=N,,0,,0 
氧 原 子 与 具 
20,(b') 2 20, (X°) + hv(762.0 nm) 


: O(S!) 2 O(D) + hv(557.7 nm) 


S t (BT): 0-0, 一 20, (b) 


NE > )+e o Na T]o«e 


辐射 率 和 临 边 柱 辐 射 强度 随 高 度 的 变化 . 图 中 , Bü 
电子 密度 的 降低 ， 气 瘤 的 体 辐 射 率 和 辐射 强度 也 随 
之 减 小 , 且 减 小 的 幅度 均 大 于 电子 密度 的 减 小 幅度 . 
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图 7 电子 (a) 和 中 性 气体 (b) 密 度 夜 间 和 白天 对 比 
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图 8 气 辉 临 边 柱 辐射 强度 和 体 辐射 率 随 电子 密度 的 变化 
(a) 电子 密度 廓 线 ; (b) 气 辉 体 辐射 率 廊 线 ; (c) 气 辉 临 边 柱 辐射 强 


E. 1 代表 原始 值 , 1/2 代表 原始 值 的 0.5 倍 , 1/3 约 为 原始 值 的 
0.33 fi 


3.2 ” 气 辉 辆 射 随 太 阳 话 动 和 地 磁 活 动 的 变化 

气 辉 辐射 来 自 于 大 气 中 电子 与 其 他 粒子 的 碰撞 、 
气体 的 离 解 复合 以 及 光 致 电离 等 作用 ,因此 太阳 活 
动 和 地 磁 活 动 对 气 辉 辐射 会 有 较 大 的 影响 . AURICH 
型 采用 太阳 射电 常数 Fio7 和 地 磁 Ap 指 数 分 别 代表 太 
阳 和 地 磁 活 动 的 强烈 程度 . 本 文选 用 两 种 应 用 最 广 
泛 的 630.0 和 557.7 nm 气 辉 进行 太阳 活动 和 地 磁 活 动 
的 影响 分 析 ， 图 9 为 它们 临 边 辐射 强度 随 太 阳 活 动 和 
地 磁 活 动 的 变化 . 

随 着 太阳 活动 的 加 剧 , 气 辉 辐射 具有 不 同 程 度 
的 增加 ， 这 是 由 于 太阳 辐射 的 增强 可 直接 或 者 间接 
使 大 气 辐 射出 更 多 的 气 辉 . 而 地 磁 活 动 的 变化 主要 
体现 在 对 气 溜 辐射 的 区 域 扰动 和 变化 上 . 当地 磁 活 
动 增强 ， 大 气 扰动 变 得 剧烈 ，630.0 nm 气 辉 的 辐射 会 
随 之 增强 , 但 其 辐射 最 主要 的 是 受 内 重力 波 等 的 影 
啊 会 有 不 同 程度 的 扰动 变化 , 这 种 变化 一 般 与 地 破 
活动 呈正 比 . 而 557.7 nm 气 辉 辐射 和 观测 的 地 理 位 
置 、 时 间 及 当时 的 粒子 密度 等 有 很 大 关系 ,变化 现象 
则 比较 复杂 ， 既 存在 正 相 关 的 情况 ,也 存在 负 相 关 的 
变化 (Leonovich 等 , 2011). 

此 外 ,基于 地 其 的 观测 结果 ,利用 模型 对 地 磁 异 
党 期 间 的 气 溜 辐射 进行 了 模拟 . 图 10 为 2007 年 1 月 16 
日 俄国 东西 伯 利 亚 地 区 在 小 型 磁 暴 期 间 观 测 到 的 气 
辉 变 化 (Leonovich 等 , 2011) 与 模型 的 比较 结果 . 图 中 ， 
实测 值 为 归 一 化 垂直 柱 辐 射 强 度 ， 而 模拟 值 为 峰值 
切 点 高 度 处 归 一 化 临 边 柱 辐射 强度 . 由 图 可 得 , COUR 
辐射 的 模拟 结果 变化 趋势 与 观测 值 变化 一 致 . 


AURIC v1.2 只 能 进行 某 一 地 理 位 置 的 单 点 气 辉 
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图 9 气 辉 辐射 强度 随 太 阳 活 动 和 地 磁 活动 的 变化 
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图 10 2007 年 1 月 26 日 俄国 东西 伯 利 亚 地 区 (52°N, 103^E) 
磁 暴 期 间 的 气 辉 辐射 实测 变化 与 模拟 结果 变化 的 比较 


模拟 分 析 ， 且 只 能 适用 于 2000 年 以 前 的 计算 , 不 适用 
于 目前 大 量 数 据 的 气 辉 模拟 及 星 上 大 气 成 分 反 演 . 
本 文 利用 Matlab 对 其 地 磁 参 数 模块 GEOPARM、 电 离 
层 IONOS 模 块 及 大 气 ATMOS 模 块 进行 了 更 新 改进 ， 
得 到 AURIC-2012. 其 中 , AEOPARM 模 块 的 改进 包括 
地 磁 经 纬度 、 破 倾角 、 太 阳 天 项 角 、 太 阳 当 地 时 间 、 
太阳 指数 Fio7 及 地 磁 指 数 Ap 等 的 计算 ; 电子 密度 的 
更 新 是 利用 国际 参考 电离 层 IRI-2012 计 算 值 代替 
AURIC v1.2 的 IONOS 模 块 计算 值 ; 大 气温 度 及 中 性 
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气体 密度 廊 线 的 更 新 则 是 利用 DTM-2013 模 型 的 计算 
IHR AURIC v1.2 的 ATMOS 模 块 计算 值 . 将 更 新 模 
块 与 AURIC 原 有 模块 相 结 合 , 不 仅 可 以 进行 气 辉 辐 
射 模拟 批量 计算 , 还 可 以 参与 星 上 大 气 成 分 的 批量 
反 演 ， 如 O/Ns 和 电子 密度 . 

将 AURIC-2012 的 临 边 辐射 强度 模拟 值 和 体 辐 射 
率 计 算 值 分 别 与 GUVI 柱 辐射 强度 测量 值 和 TIDI 体 辐 
射 率 测量 值 进行 比较 ,辐射 强度 峰值 模拟 平均 相对 
误差 为 18.6%， 体 辐射 率 峰 值 模拟 平均 相对 误差 为 
15.45%. 最 后 ,利用 AURIC-2012 进 行 多 年 份 春分 时 
的 昼 气 辉 和 夜 气 辉 临 边 辐 射 强度 随 纬 度 和 高 度 的 二 
维 模拟 ， 并 基于 模拟 结果 对 昼 气 辉 和 夜 气 辉 的 辐射 
强度 及 分 布 规律 进行 了 比较 分 析 . 由 模拟 分 析 可 得 ， 
春分 点 附近 氧 原子 、 氧 气 分 子 及 氮气 分 子 气 辉 辐射 主 
要 位 于 中 低 纬度 地 区 <500 km 高 度 范围 内 ,， 随 着 纬度 
的 增加 ,辐射 逐渐 减弱 ， 而 随 着 高 度 的 增加 辐射 一 般 
先 增 大 后 减 小 . 此外, 不 同 谱 线 气 辉 的 辐射 强度 差异 
较 大 ,同一 谱 线 气 辉 白天 和 夜间 的 辐射 差异 也 较 大 ， 
这 些 现象 是 由 气 辉 本 身 的 辐射 机 制 、 大 气 粒 子 的 密 
度 、 太 阳 活 动 和 地 磁 活 动 强度 等 共同 作用 的 结果 . 

气 辉 的 批量 辐射 模拟 分 析 , 不 仪 为 星 载 测量 提 
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供 物 理 设计 理论 指导 ,也 为 星 载 大 气 成 分 探测 提供 ” 体 成 分 进行 密度 更 新 计算 ， 同时 进行 AURIC-2012 
hd B. 但 由 于 中 性 气体 模型 代码 的 限制 ， 辐射 物理 计算 模块 的 珍 代 ， 并 在 此 基础 上 对 茶 些 气 
AURIC-2012 只 对 五 种 主要 中 性 气体 成 分 进行 了 密度 ” 瘤 的 辐射 机 制 进 行 修正 和 参数 更 新 ， 以 更 准确 的 进 
更 新 ， 因 此 在 今后 的 工作 中 需 进 一 步 对 其 他 中 性 气 。” 行 更 多 气 辉 谱 线 的 辐 册 模拟 和 星 上 数据 反 演 研 究 . 
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